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We propose a variational perturbation method based on the observation that eigenvalues of each 
parity sector of both the anharmonic and double-well oscillators are approximately equi-distanced. 
The generalized deformed algebra satisfied by the invariant operators of the systems provides well 
defined Ifilbert spaces to both of the oscillators. There appears a natural expansion parameter 
defined by the ratios of three distance scales of the trial wavefunctions. The energies of the ground 
state and the first order excited state, in the 0*'' order variational approximation, are obtained with 
errors < 10~^% for vast range of the coupling strength for both oscillators. An iterative formula 

, is presented which perturbatively generates higher order corrections from the lower order invariant 

' operators and the first order correction is explicitly given. 
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There are two well known problems in studying the anharmonic and the double-well oscillators. The first is that a 
naive perturbation theory of anharmonic oscillator leads to a divergent perturbation series even for an infinitesimal 
I value of the coupling strength due to the eventual dominance of the perturbative correction over the un-perturbed 
' contribution for large amplitude of oscillation [l| . To remedy this difficulty several generalized perturbation methods 
T-H , have been developed 0, 0| ■ There have been many recent developments which provide well defined perturbative 
CN ■ iterative solutions for the eigenvalues and eigenfunctions 0, 13, S B [3l ■ interesting differential equation method, in 
, parallel with Mathieu equation, was presented by Liang and Miiller-Kirsten Q . A through analysis for this asymptotic 
perturbation and large order behaviors was given in Ref. . 

Another difficulty is that a double-well oscillator has completely different algebraic structure and ground state 
wavefunction from those of anharmonic oscillator. Therefore, most perturbative approaches treat the two oscillators 
separately 0, 0]. The transition from a single- well system to a double- well system of a time-dependent oscillator 
was calculated through a nonperturbative method such as the WKB approximation [Til , [l^ or by^ accommodating 
^ ' contributions from many excited states in addition to the ground states of symmetric oscillator 13]. Even though 
these methods provide good approximations to the ground state and excited state energies, its generalization to 
field theory of infinite degrees of freedom or to time-dependent systems with symmetry-breaking-like phenomena 
encounters difficulties. Therefore, it is important to establish a consistent approximation method that can treat both 
of the oscillators in the same way. This subject has not been extensively studied compared to the problem of divergent 
Ci ] perturbation. 

The purpose of the present paper is to present a possible resolution to the second problem using a variational 
perturbation method. We present a unified treatment of the anharmonic and the double-well oscillators in the point 
of view of a generalized deformed oscillator. This provides an algebraic method in which the two oscillators are 
dealt with on an equal footing. We also present a natural expansion parameter defined by the ratios of three length 
scales appearing in the trial wavefunction. Its numerical value is smaller than 1/2 for all physical parameter range 
of the system for the double-Gaussian trial wavefunctions used in this paper. To do these, we solve the Liouville-von 
Neumann (LvN) equation perturbatively. 

The generalized deformed oscillator 14 1 was developed in the approach of quantum group theory [isl . [l6j . The 
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algebra is generated by the operators {1, a, a^, N} and the structure function ^(x), satisfying the relations 

[a,N] = a, [a\N] = -a\ (1) 
a^a = $(iV) = [N], aat = $(iV + 1) = [A^ + 1], 
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where N is the number operator and ^{x) is a positive analytic function with the condition 

$(0) = 0. 

This ensures the vacuum state |0) to be defined as the eigenstate of N with zero eigenvalue and satisfies 

a\0) = 0. 

It can be proved that the generalized deformed algebra provides a Hilbert space of eigenvectors |0), |1), ■ ■ ■ ,\n), - ■ ■ 
of the number operator N with (n|m) = 5mn- These eigenvectors are generated by the formula, 



V L'-J- k=l k=l 

The operators and a are the creation and the annihilation operators of this deformed oscillator: 



a\n) ^ ^/\n\\n - 1) , a^n) = [n + l]\n + 1) . 

q" - g"" 

The harmonic oscillator is obtained for $(n) = n and the deformed oscillator for <i>(n) = — j- with the 

Hamiltonian, H = -^{aa^ + a^a). 

The Hamiltonian describing both the anharmonic oscillator and the double-well oscillator is given by 

with being a positive (anharmonic) or a negative (double-well) real number. For A — > with u'^ > 0, the system 
becomes a harmonic oscillator (harmonic oscillator limit) and for < with A ^ 0+ the system describes a set of 
two separated wells of harmonic oscillators (infinitely separated double- well limit). 
The LvN equation [13], 

dA{p,x,t) dA ^ 

is equivalent to the Schrodinger equation and physical information can be obtained from the invariant operator A. 
Based on this operator, a complete Hilbert space of the oscillator can be constructed. The equation ^ was used 
as a starting point in developing many approximation methods and of physical applications (isl . [l9j . This method 
provides a powerful basis in constructing the variational perturbation methods since the zeroth order solution of ([3]) 
can be made to be a variational approximation of the system. In the present paper, we will utilize this aspect to 
develop a new variational perturbative approximation with two variational parameters for time- independent oscillators. 
Generalization to time-dependent systems is straightforward since the LvN equation is known to be one of the best 
ways to treat the time-dependent system [itI . [isj . 

Standard perturbation theory based on harmonic oscillators assumes that the coupling term Ax'* is smaller than the 

other term: (Ax*) <C (uPSy^). This leads to a 'localization condition' of the expectation value, iSp') ~ — <C w^/(3A), 

with respect to the ground state wavefunction. This inequality restricts the applicability of the standard perturbative 
method to small coupling. To overcome this limitation, the Gaussian approximation method [13] and its generalizations 
have been developed p3. [25j. The Gaussian approximation starts from defining a new frequency scale Vtc, which is 
variationally determined to satisfy 

= ^2 . 3A 



3A 

The localization condition now takes the form ^ 1, which is satisfied for uo"^ > 0. The characteristic length scale 
and momentum scale corresponding to this frequency are 
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The Gaussian wave-function is annihilated by the annihilation operator, 

ip X 
2pa 2xg 

which satisfies a simple commutation relation [acOc] ~ ^- Based on this commutation relation, a complete Hilbert 
space can be constructed by consecutive operations of the creation operator to the ground state |0)g defined by 
aG|0)G = 0. 

There have also been some attempts to apply the Gaussian approximation to the double- well systems, in which one 
wants to find an annihilation operator of the form 

. _ _ip_ x-xq 

2pG 2XG 

where xq and are the parameters to be determined by variation 0, 0| • In spite of the success of this approach 
in finding an approximate ground state energy of double-well oscillator, it is to be noted that the approach does not 
reflect the symmetry of the system in the sense that the state annihilated by dx^ does not have the symmetry of 
the potential, the space inversion symmetry. An approximation method which keeps this symmetry is the Double 
Gaussian Approximation (DGA) [8l.T2ll.[23j. in which the sum of two Gaussian wave-packets (or similar ones) are used 
to approximate the ground state wave-function of a double-well oscillator. A major deficiency of the DGA is that no 
Hilbert space based on the approximate ground state has been found, which is one of the issues to be dealt with in 
the present paper. 

We develop a method in which a double-well oscillator is dealt with in the same way as a single-well anharmonic 
oscillator. Before introducing the method, we point out some interesting features of double-well and anharmonic 
oscillator systems. The first is that the eigenvalues of a double-well oscillators are not approximately equi-distanced 
contrary to the case of an anharmonic oscillator. For example, the difference between the first excited state energy 
and the ground state energy is considerably different from the difference between the second excited state energy and 
the first excited state energy. This implies that a single set of creation and annihilation operators (a, a^) cannot 
connect all the eigenstates with reasonable accuracy. Considering the parity even states and odd states separately, 
however, one finds that the energy eigenvalues are almost even-spaced in each parity sector [For example, see chapter 
5 of Ref. [13] for exact solution of double-well oscillator with the potential ^(fclxj — a)^]. This suggests us to consider 
operators connecting eigenstates by two steps, connecting the states among the same parity eigenstates. We call the 
ground state |0) and the 1** excited state |1) as the 'ground' states for each parity sector. Thus, the f** excited state 
(|0o) = is the ground state of odd parity sector and the state jO^;) = |0) is the ground state of even parity sector. 

Another interesting feature of double- well oscillator is that there exists a set of operators {Aq, Aq) and even and 
odd ground states {\Oe/o))j which consistently describe both the harmonic oscillator limit and the infinitely separated 
double-well limit. Consider the following states of even and odd parities, 

|. . 1 1+) + !-) |. . 1 l+>-|-> 



Ne V2 ' ' ' No V2 ' 

where the states |±) stand for the ground states of a harmonic oscillator centered around x — ±xo, annihilated by 

ip x^xq 
a± = h — , satisfymg 

2pG 2XG 

a±|±)=0, a±|T)-T|T) — , 

XG 

and the normalization constants are Nq = f — (+|— ) and Nj^ = f + (+|— ). The time-independent wavefunctions for 
states |±) are given by the Gaussian wave-packet. 



vl/±(x)^(x|±)= (^j e-^^^^^^'^^l^) e . (4) 



The overlap of the two states |±) used in defining the normalization constant is 

(+|-> = e^T^^^. 

Evidently, these even and odd ground states are orthogonal to each other: 



o(0£|0o)o = 0. 
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We now show that the states \Qe/o)o ^nd an operator (AJ ~ ^^a^a^), where x ~ ^JJx^, correctly describe the 

two hmiting Hilbert spaces, the harmonic oscillator limit [uP > 0, A — > for xq — > 0) and the infinitely separated 
double- well limit {lo^ < 0, A ^ for xq — *■ oo) of the system described by the Hamiltonian 

Consider the harmonic oscillator limit first. The ground state wavefunction of a harmonic oscillator is a Gaussian. 
Since the sum of two Gaussian wave packets of the same frequency is a Gaussian, the ground state wavefunction is 
correctly reproduced in the xq — s- limit. Let us examine the 1^* excited state in the xq — > limit. The odd parity 
ground state becomes an exact first excited state wavefunction of a harmonic oscillator in the xq = limit: 

^hm (x|0o)o - (^) /^-exp . 

Therefore \Qe/o) gives the exact ground state and the first excited state in the harmonic oscillator limit. Repeated 
actions of the operator ^ a^aL ~ a^^, with xq = x, on the two states 1 0^/0)0 reproduce the complete Hilbert 
space of the harmonic oscillator. Both of the states \Qe/o)o ^re annihilated by ^ a^a_ = d^. In this sense, the set 
of operators {Aq, Aq) and the even and odd ground states \0e/o) correctly reproduce the harmonic oscillator system 
in the xq ^ limit and Aq and Aq play the role of an annihilation operator and a creation operator which raises and 
lowers the states by two steps. 

We next consider the infinitely separated double-well limit by setting xq — *■ 00. Now, the two states |±) are 
completely separated from each other and do not overlap. Therefore, the energy of the two states \Oe/o)o 
degenerated and the operator Aq acts separately on each states |±) and raises each states |±) to the l'** excited 
states. For example, for |+) state, 

-^alal\+) = —a\ [^1 + — J 1+) = y«+l+> + O(^) ^ |1+), 

where in the infinitely separated double- well limit, xq/x ^ 1 and xq/x 0. Similar calculations will be applicable 
to the state |— ). In addition, Aq annihilates both of the states |±), simultaneously. In this sense, the operators Aq 
and Aq play the role of an annihilation operator and a creation operator and correctly produce the Hilbert space of 
the infinitely separated double- well oscillator in the Xq ~^ 00 limit. 

This inspection of the two limits xq ^ 0, 00 implies that, by using the even power operators ^0 a-nd Aq, one may 
describe both the anharmonic and the double-well oscillators simultaneously. 

We construct a 0*'' order invariant annihilation operator A of the exact annihilation operator A of the system ([2]) 
in Sec. II by solving the LvN equation to the lowest order for variational approximation. We construct the creation 
{A^) and annihilation (A) operators as even-power series in p and x in such a way that they raise or lower the energy 
eigenstates by two steps. For example, the creation operator raises the ground state to the 2"'' excited state and the 
1** excited state to the 3'''^ excited state, etc. The annihilation operator A must annihilate both of the ground state 
jO^;) = |0) and the first excited state |0o) = and their wavefunctions are an even function and an odd function 
of X, respectively. In the harmonic oscillator limit, xq — > 0, the annihilation operator must become the square of 
the annihilation operator of the harmonic oscillator. In the infinitely separated double-well limit, xq — > 00, the 
annihilation operator A must reproduce the annihilation operator of a harmonic oscillator which is centered at each 
bottom of the double-well potential. An algebra satisfied by the operators are obtained and then used to construct 
the Hilbert spaces for the oscillators. This operator A is used in Sec. Ill to construct the variational approximations 
of the anharmonic and the double- well oscillators. It is shown to provide a good approximation for the energies of the 
oscillator eigenstates with errors smaller than 10~^% for most range of parameters of both the double-well oscillator 
and the anharmonic oscillator. We then calculate the first order correction to the invariant operators {A,A^), and a 
systematic method to construct higher order invariant operators in perturbative series is developed in Sec. IV. We 
conclude with some discussions on our method in the last section. 



II. ZEROTH ORDER SOLUTION OF THE LIOUVILLE-VON NEUMANN EQUATION 

The variational perturbation method starts from variationally identifying the 0*'' order approximate ground state 
and an annihilation operator which annihilates the ground state. For example, the Gaussian approximation uses the 
Gaussian wave-packet, Ne~^'^^ as the zeroth order ground state, with VI a variational parameter. In this section, 
we try to find a 0*'' order ground state and an annihilation operator, which describe both the double-well oscillator 
and the anharmonic oscillator, consistently. 
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A. Localization condition 



There have been attempts to find the solution a of the LvN equation ([3]) for anharmonic oscillator as a Taylor-like 
series of p and a; [l^ . This process reveals that the anharmonic oscillator ^ with lo^ > Q has the structure of a 
g— deformed oscillator to the first order in the variational perturbation series and that of a generalized deformed 
oscillator to the next order [2^. This observation is based on the fact that the low lying eigenstates of anharmonic 
oscillator are localized around a: = and p = 0, in the sense that, 

(i>-0, (i'>-^, (?5)^0, {f)-^. (5) 

In the Gaussian approximation, this condition is implicitly assumed. But for the double- well systems the conditions ([5]) 
are not satisfied. 

The probability distributions of low-lying eigenstates of a double-well oscillator would be bi-centered around x — 
iba^o (2:0 > 0) with {p) — 0. We thus develop a generalized series expansion in terms of bi-centered polynomials. For 
example, a non-singular function /(x) can be written as 



/(-) = E 



(6) 



Identifying for n = 0, 1, 2 • • • provides an approximation of f{x) around x ~ ±a;o in a series form. If one wants to 
describe physics around x — ±xo, this approximation can be used effectively. 

We assume the energy scale of each packet at x — ±xo to be hTl. Using dimensional analysis, we have the 
'localization condition', 

(x)^0, {ix'-xir)^{2xr^, (p)~0, ip')^^, (7) 



where x = ^J{SF) . In this condition, three different length scales appear: the position of each bi-centered packets 
(xo), the root-mean-square expectation value of the position operator (i), and the length scale corresponding to the 
energy scale {xq_ = ^Jh/{2^)) . The ratios of the three scales provide interesting dimensionless non- negative quantities 
which are not larger than 1, 



The inequalities can be easily proven for the wave-packet given by the sum of two equal packets: 

1 

''P±{x) = [(l){x — xq) ± (l){x + xq)]; N^1± dx'4>{x' — xo)4>{x' + xq)- 



/2N 

With respect to this packet, we have the relation 

2 

-2 2 I Xq 



X 



N' 



where we have assumed the wave-function being real for simplicity. For ip- case, Eq. ([8]) is evident. For ip+ case, the 
value e = Xq/x'^ will be maximized at xq 00, where e = 1. Therefore, the inequality in Eq. ([5]) is valid. 

We consider an explicit example of a bi-centered packet given by the sum of two Gaussian packets centered at 
x = ±a;o, respectively. Then we have the inequality x > xq and x > xn always, where the equalities hold for xq = 00 
and xo — 0, respectively. Thus the ratios e and e are always not larger than 1 and they take the values e = l(e = 0) 
for xo — Q and e = 0(£ = 1) for xq ^ 00. Since the multiplication of these parameters, 

5 = ee, (9) 

is always smaller than 1, it can be used as a good perturbation parameter for systems with bi-centered wave packets 
as approximate eigenstates. Note that these parameters are not fixed, but depend on the nature of wavefunctions. 
Thus we will use these parameters as variational parameters. 

For calculational convenience we introduce dimensionless operators tt and y for momentum and position operators 
defined by 

^ ' y-^- (10) 
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FIG. 1: Schematic plot of the values of e, e, and 5 as functions of the distance between the central positions of the two Gaussian 
packets. This is the explicit zeroth order result for even-parity eigenstates. The horizontal distance in the figure is rescaled by 
the tangent function, j/o = tanz. Note that the explicit numerical value of 5 is smaller than 1/2. 



Note also that yo — xq ^JVl/H relates the two parameters e and e by 



The commutation relation between y and tt is 



and the localization condition ([7]) becomes 



2^yl = e. 



The expectation value of is explicitly given by 



1 



For later use, we introduce some operators of 0(e") and 0(e°), 

T^ = V2fr, T^^V2ey, T^ = 2it\ ^ e{f - yl), TI = e{T:,y}+, = e\f - yl)\ 
which will be useful in describing the bi-centered systems. 



(11) 
(12) 

(13) 
(14) 
(15) 



B. Hamiltonian for variational approximation 

We now consider a quartic oscillator, which will be used as a basis for the variational approximation, described by 
the Hamiltonian: 



Hn = ——H, 



Vo; 



H„ 



1 



* rinO , ^^rp2 , ■'^^ rp4 



(16) 



where Vq is the potential value at y = yo and the subscript in Hq denotes that this Hamiltonian will be used as a 0*'' 
order form for the variational perturbative expansion of double-well and anharmonic oscillators. The constant eo is 
the value of e in the harmonic oscillator limit, y^ = 0. Explicitly, eo = 1 for \Qe)q a-nd eo — 1/V3 for |0o)o- Note that 
this is one of the most general form of quartic Hamiltonian if $7 and yo are chosen arbitrarily. Note also that only one 
of the three quantities e, e, and yo is an independent variable because of Eqs. (O, (fTTj) . and p4|) . Since we are to use 
the ground states of (|16p as trial wavefunctions for the variational approximation for the Hamiltonian ([2]), we have 
two variational parameters, 17 and yg. 
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The coefRcients {v2, A) in the Hamiltonian ([T6)) can be written in terms of the above mentioned parameters by 
considering the two hmiting cases, the harmonic osciUator (j/q — 0) and the infinitely-separated (j/o oo) double-well 
limits. 

In the harmonic oscillator limit, the ground state is given by a Gaussian packet. The exact frequency is given by 
= uj and there is no quartic potential in the Hamiltonian. Therefore, the Hamiltonian (jl6|) becomes that of the 
harmonic oscillator with correct frequency lo if 

harmonic oscillator limit (e ^ 0, e — > en) : I I'' 

^ ' [ A(j/o ^ 0) = finite number. 

On the other hand, for the ground state of the infinitely separated double-well limit (j/o ~* oo)j we have e = and 
e — \ and the Tq term disappears since e vanishes. Each packet at x = ±xo evolves as if it is a free harmonic oscillator. 
The frequency of each packets is determined by the second derivative of the potential with respect to x at the bottom 
of the potential. This determines fi^ — 2Axq = 2\x^ . Therefore, the coefficient A of Tq term must be 1 in this limit: 



infinitely separated double- well limit (e 1, e ^ 0) : 



/ ^2(2/0 ~* 00) = finite number, 
\ A(yo ^ 00) = 1. 



With these observations, it is clear that the constants V2 and A are of 0{5^). The values V2 and A become important 
at the harmonic oscillator limit and the infinitely separated double-well limit, respectively. For both limits, therefore, 
we may set 

^'2 - 1 = A, (17) 



for the zeroth order approximate Hamiltonian Hq. With the values of V2 and A of Eq. (|T7l) . the Hamiltonian (|T6|) still 
describes the most general anharmonic and double-well oscillators since O and uq are free parameters. 




C. Zeroth order results 

As we have discussed in the previous section, we try to find an annihilation operator A and a creation operator 
to 0{5^) in this subsection. Since we want to find the operators as even functions of the bi-centered operators we try 
the simplest ansatz for the zeroth order operator A: 

A - u{t)T^ + v{t)Tl +w{t)T^. 

dAip, X, t) dA 

To the zeroth order, one can easily find that the LvN equation, ih ^ — ^ — ~ ifi-^ — h [A, Hq] = 0, becomes a matrix 

at at 

differential equation of the form 

-^ = -MAit), Ait)- 

where Af is a (3 x 3) matrix determined by the commutator of A and Hm of Eq. (|16p . For the time- independent 
case, the differential equation can be solved by a simple matrix eigenvalue equation by setting A(<) = e"'^*A(t = 0). 
Explicit calculation of the commutator [A, Hm] gives the matrix M: 

2eVe2 + ^2 ' 
M = -i \ -e 
-1 2e/el 

The eigenvalues of M are given by the solution ( of 
where 



c2 



(18) 
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We have three different C eigenvalues, (0,±a). The ( = solution is, in fact, the Hamiltonian itself. Since what we 
want to find in this section is the creation and the annihilation operators, we do not try to analyze it here. 

The other two solutions with eigenvalues C = ±a denote, in fact, a set of operators A and A'' . If we demand A 
being a negative frequency mode, we have 



A = Doe a; a = —iv 1q + H 

A' = Doe a, a' — iiy 1q + li 

a 

where bo is a normalization constant to be determined later and v is 

^.2 2eVeg + £^ _ 
a 



(19) 



(20) 



These are the zeroth order annihilation and creation operators we were looking for. Thus the operator A satisfies 

dA{t,p{t),x{t)) 



dt 



0, 



where ~ denotes "the same up to 0((5'')" and A and .4^ are the zeroth order invariant operators. The explicit algebra 
satisfied by the operators A and A^ will be shown in the next subsection. 

D. Algebraic properties of the zeroth order invariant operators 



In the previous subsection, we obtained the annihilation and the creation operators which satisfy the zeroth order 
LvN equation. We expect that the zeroth order solution of the LvN equation has the structure of a generalized 
deformed oscillator as we have pointed out in the introduction. In this subsection, we show that there exist well 
defined algebraic properties of the operators A and A^ . 

Using Eqs. and IT^ . we find the commutator of the operators A and A^: 



= V e 



ae 



a 
a 



OiS). 



(21) 



The anticommutator of the operators A and A'^ is 



46§ 



o 2 

~2~ 



2e2 fl 



-T" 



(22) 



0{5). 



The algebraic structure of the deformed oscillator will be apparent by introducing a structure function $o(-^o) of 
the form 



$o(iVo) = Abli^E+No + Ae'Ni 
, and (ESI) with Eq. ([1|), we have 



(23) 



By comparing Eqs. ([21 



The constant £q is determined by the condition $0(0) = 0. With (pij) and 

2 



the condition reads 



ae \ 1 

a 



2e' 



So 



1 



(24) 
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This equation determines the energy eigenvalues in both the harmonic oscillator and the infinitely separated double- 
well limits. For the infinitely separated double- well limit (?/o — ^ cxd, e — > 0), we have £q = a and for the harmonic 
oscillator limit (j/o ^ 0, e ^ 0), we have two different solution 8o = 3a/2, a/2. These two eigenvalues correspond to 
the energy eigenvalues of the odd-parity and the even-parity ground states. The rescaled ground state energy Eq in 
Table 1, which is determined by Eq. ((59)) below, satisfies this condition l|24[). 
The structure function (j23[) and the number operator Nq, 

No = (25) 
satisiy the conditions for the generalized deformed oscillator, Eq. |T]), since 

iU~$o(^o), AA^ ~MNo + l). (26) 
The Hamiltonian can be written in terms of the number operator Nq as 

H ^^{2aNo + £o) + Vixo). (27) 

nn 

Therefore, the ground state energy is given by ~^^o + ^(a^o) and the energy difference between two nearby states in 
each parity sector is ahfl. 

The yet undetermined quantity bo in Eq. (fT9|) is fixed by imposing [.4,.4^]|Ob/o)o — \OE/o)a- This gives 

Note also that due to the LvN equation ITBl) . and the definition of No we have the desired commutation 

relation: 

[A,No]~A, [A\No]~-Al (29) 
With these algebra and Eq. ([23]) . the ground state \Oe/o)o can be defined as an eigenstate of No'. 

A\Oe/o)o - 0. (30) 
From Eq. (|26p and the condition $o(0) = 0, the number operator annihilates the ground states: 

No\Oe/o)o = |OB/o)o<fo '(0) = 0. (31) 
The n*'' excited states are generated by the formula: 

\nE/o)o = ^^\Oe/o)o, (32) 

where 

n n 

[n]ol^l[[k]o = l['Po{k). (33) 

k=l k=l 

The states \nE)o and \no)o are the {2rif^ and the (2n -I- 1)"' excited states with even and odd parities of the 
anharmonic and the double-well oscillator. Due to Eqs. (pS)) and the states have the eigenvalues tie/o for iVo: 

^o|'^£;/o)o - nE/o\nE/o)o- (34) 
Moreover, these eigenstates are orthonormal to each other to the zeroth order: 

o{nE/o\mE/o)o - Snm- (35) 

Eqs. (|23|) . ([25]) . (|26|). ((29|) . and ([30]) constitute the conditions for a generalized deformed oscillator which possesses a 
Fock space of eigenvectors \Qe/o)oi |1_b/o)oi ■ ■ ■ j |'T'E/o)oj ■ ■ ' of the number operator A^o- Further specification of A 
to higher order is possible perturbatively, which is the purpose of Sec. IV. 
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III. VARIATIONAL APPROXIMATION 



In the previous section we have obtained a zeroth order invariant creation and annihilation operators which satisfies 
the LvN equation for the zeroth order Hamiltonian (jl6p which depends on two arbitrary parameters, fl and j/o- By 
using the ground state of the zeroth order Hamiltonian (jl6p as a trial wavefunction we develop, in this section, a 
variational approximation of the anharmonic and the double-well oscillators. We then compare our result with those 
of other approximation methods. 

In coordinate space y, the eigenstate vE'(y), which satisfies A'i'{y) = 0, is given by the linear combination of 
g-"y Hm{f3y) and the confluent hypergeometric function iFi, which give difficulty in analytic manipulations. There- 
fore, we propose to use an operator clq which is equivalent to A up to the zeroth order but have simpler ground 
state: 

isS n A ^ 

where 

a± = ^^ + ^(yTyo) (36) 

are the annihilation operators of harmonic oscillators centered around y = ±yo- The new operator is explicitly written 
as 

ao = -^^''Tl+Tl + ^T^. (37) 



The operators tt and y are related to a± by 



2v 



^{d±-ai), (38) 



vyo = a+ - a^, 



Thus A is equivalent to the zeroth order to the product of the two annihilation operators of harmonic oscillators 
centered around y = ±yo- Note that the annihilation operators a± and the creation operators a\. satisfy the standard 
commutation relation 

[a±,al] = l, [a+,a_]=0. (39) 
Consider the states |±) annihilated by d±, 

a±|±)=0. (40) 

The wave- functions for states |±) are 

.1/2 



*±(y) ^ {y\±) = — e-"- (^-^°) '\ 



which determine the overlap 



2 2 



Since both of the two states |±) are annihilated by do, the approximate ground states of even or odd parity, \Qe/o)i 
which are annihilated by do, are given by the linear combinations, 

i„ ^ 1 1+) + !-) |. ^ 1 1+)-!-) 

The normalization constants are 

Ml ^ I + {+\-) = 1 + er^'yo ^ Ml^l-{+\-)^l-e-''"y». (42) 
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The ground states (|4T|) satisfy 



o("-E/o|0£;/o)o = 0, for 717^0, 



where the excited states \nE/o)o a-re generated by the successive action of a\ on the ground states. We thus have 

constructed the set of operators (ao,aJ) and the even and odd ground states (|0£;)o, |0o)), which properly describe 
both the anharmonic and double-well oscillator simultaneously to the zeroth order in 5. 

The expectation values of operators with respect to the states \Qe/o)o can be obtained algebraically by using 
Eqs. and For example, the expectation value of is 



o{(^E/om\OE/o)^ = -^{{+\±{-\){a+-a\Y{\+)±\-)) = u'TE,o, 
where upper (lower) sign refers to the even (odd) parity case and 



e-'vi ± 1 ' 



(43) 



(44) 



Starting from 1, Te decreases to 0.45 at i/i/q = 1.1 and then increases again to 1 as vyo increases. On the other 
hand, Tq continuously decreases from 3 to 1 as vy^ increases. The functions Ce/o are 



Csivyo) 
Co{fyo) 



-)/Ml^{e'''yi + iy\ 



iviE 



2 2 

Vo 



As shown in this equation, the expectation values for the even and the odd ground states take similar form except 
for the indices in Ce/o- From now on, C will denote both of Ce/o if not specified. This convention is applied to all 
other variables such as e, v, M, Eq, T, and R defined below. The expectation value of Tg is 



o(0 



E/O 



^0 |Ob/o)o 



The formal identity x'^ = {{x^ — Xq)) + Xg leads to the following relation among parameters. 



AP 



Ho? 



(45) 



(46) 



Since the normalization constant M is a function of vy^, this equation determines ej/g in terms of vy^ or vice versa. 
This equation has interesting implications on the limiting cases. In the limit t/o ~* 0, Eq. (j46p determines = 1 and 
Co = l/\/3. On the other hand, in the limit j/o — *■ oo, we have e ~ 1/(a/22/o) irrespective of the parity of the states. 
The equation (|46p can be rewritten for e as a function of i/g = vya'- 



1 + 



(47) 



In addition, e can be written as a function of j/g, from Eqs. (|18p . (PO)) . and (|47p . as 



£(yo) 

giv) 



V3 
V2 



3^0 + 1 , 1/3 



(48) 



1 



+ iViyl^l + llyt~yl 



The explicit values of e, e, and 5 — ee are plotted in Fig. 1. 

We also calculate the expectation value of the quartic part of the Hamiltonian: 



o(0 



e/o\Tq\Qe/o)o 



3 + 8iy2j;2 



4P* 



3 + Sy'^yl ± 



+ v^yl (v^yl - 3) C 
^>^hl{v^yl-i) 



(49) 



e^'^vo ± 1 
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A. Effective Potential 



p2 ^2 

As in the Gaussian approximation, in which the ground state oi H = ~2^~2'^ chosen to be the trial wavefunction 

for the Hamiltonian we use the states ([1T|) . which are the 0{5^) even and odd ground states of the Hamiltonian (|16p . 
as trial wavefunctions for the variational approximation for the Hamiltonian Therefore, the yo ^ Q limit of our 
approximation becomes the Gaussian approximation, and the approximation with yo — l/(V2e) limit describes a 
double Gaussian approximation 8, 21]. For this variational approximation we write the Hamiltonian ([2]) of a general 
anharmonic llP' > Q or a double- well cj^ < oscillator in terms of operators defined in Eq. (fTS]) : 



H 



— + —SP- + -x^ + —t 
2 2 4 4A 



(50) 



1, 



Ml 



1 f UJ 



o-'2+lr>2+ o2 Mo 



Xh 



rp4 

J- n 



2f]3e2 



1^0 , 1 sign(cj2) 2A^yg 



+ 



X. 



rpl 

n 



where the potential V{x{)) is 



sign(c.^)^ + 



(51) 



In this equation the theta function, 9{-~liP'), is included to make the lowest energy of the potential vanish for lu'^ < 0, 
and xq = \J h/Vl 2/0, and the dimensionless parameters and Vlj. are defined by 

Xh o - ^ 

2k7F' ^"M- 



At- 



(52) 



From Eqs. ([33]), ([jS)) . and the expectation values of the Hamiltonian (|5n|) with respect to the ground states 
|0_e) and |0o) in Eq. (gT]) become 



2{H) nrv'^T . , 2^T + 2z/2y2 

where we have used Eq. (|46p to ehminate e for t/, and 



At- 



4(f7^Z/2)2 



4Xr ' 



i?(i.yo) - 3 + 12i.2y2 + 4^4y4 _ 4^2^2(3 ^ ^2^2)^ ^ 3 ^ 



4^2^2(3 ^^2^2) 

1 ± e-^''?'o 



The function Re {Ro) monotonically increases from 3 (15) as vyo increases. Note that the rescaling. 



^yo = yo : 



(53) 



(54) 



(55) 



completely remove the z/ dependence on the expectation values since T defined by (|44p is the function only of vyo. The 
value of only affects the relation between e and yo through Eq. With these new parameters, the expectation 

value (l53l) becomes 



2{H) riTiyo) , . , 2^7M±Ma. R^- ^ 



-0.2) 



4A^ 



The variation of (f56|) with respect to £7 gives the gap equation 

- sign(u;2) ( 1 



Tim) 

which determines non-trivial unique positive solution for Q 



2yl \^ XrRiyo) 



0, 



fi(yo) = ( 3 



1/3 



3~2V3 



sign(a; 



2^0 



Di/3 



(56) 



(57) 



(58) 



D 



9Xr.R{yo) 

T{yo) 



2\/3W-sign(w2) ( 1 



T(yo) 



27 A2j?2(^„) 
4 T2(yo) ' 
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for all values of uJ^ irrespective of its sign if A > 0. We set this solution be fi(yo)- Note that the gap equation ([57|) 
takes similar form as that of the variational Gaussian approximation except for the dependence of the coefhcients on 

We thus have the effective potential as a function of j/o by inserting the solution ([55]) to Eq. 



2K//(yo) 



f^(yo)r(yo) 



:^(yo) 



h\uj\ An^yo) ' 4A, 

The functional dependence of Veff{yo) on yo is shown in Fig. [2] for a couple of values of Ar. 



(59) 



. 561 
0.5605 
. 56 
0.5595 




1 . 55 
1 . 54 
1.53 
1.52 
1.51 




yo 



0.2 



FIG. 2: Effective potential ^'^■f-f g, function oiyo for the coupling constants Xr — 2 (left) and Ar — 20 (right) for o;^ > 0. 

h\u;\ 

The energy decreases very slowly for yo <^ 1 and increases for large yo- Therefore, there exists a unique minimum of 14//(j/o) 
for positive yo- Since the graphs are nearly flat for yo <^ 1, yc increases very fast for small Ar as depicted in Fig. 3. 



We write the value yo that minimizes the effective potential (|59p as j/c 

dVeffivo) 



dyo 



0. 



(60) 



Then, the ground state energy is given by Veff{yc)- With this y^ ^{yc) is determined from Eq. and e(t/c) 

and e(yc) are determined from Eqs. (|47)l and (|48|) . With these results the annihilation operator ([37)1 is uniquely 
determined. In Appendix A, we explicitly analyze the variational equation, Eq. (|60p and the results for positive 
are shown in Fig. 3. 



1.4- 
1.2- 
1 - 
0.8- 
0.6 - 
0.4- 
0.2- 




15 



15 



20 



FIG. 3: Behavior of ground state energy ( — -rA-r—)^ 3^=' '^^'^ functions of Ar (horizontal axis) for Lo^ > 0. 

The left figure shows the energy (grey curve) and j/c (black curve) as functions of Ar. On the right shows the behavior of the 
present f2(Ar) (black curve) compared to QciXr) (grey curve), the Gaussian value. As seen in this figure, the effective frequency 
n is much larger than the Gaussian value Qq for Ar 3> 1. 

In general we have nonzero yc for positive A. We present the energy eigenvalues of the even and odd ground states 
and £ for several values of the coupling strength in Table 1. As one can see in Table 1, even for a very small coupling, 
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such as A,- ~ 2 • 10"'^, we have non-neghgible j/c value, e = yl/y'^ ~ 0.1. Since we have ?/c = in the Gaussian 
approximation, non-vanishing measures the departure from the Gaussian results. This is measured by the value e 
in Table 1. The energy eigenvalues of our zeroth order approximation is closer to the exact values than the 2nd order 
perturbative results of the approximation methods based on the Gaussian approximation for vast range of couplings. 



TABLE I: The explicit energy eigenvalues for anharmonic and double-well oscillators for several values of coupling strength. 
'Caswell' and 'NGAS' refer to the results of Ref. [23] and Q], respectively. 'Exact' implies the exact numerical result. 





oj^ > 


<Q 




0.002 


0.02 


0.2 


2 


20 


200 


200 


20 


2 


0.2 


e 


0.1054 


0.2824 


0.5170 


0.6419 


0.6751 


0.6824 


0.6865 


0.6940 


0.7296 


0.8705 


Eo/{huj) 


0.5007 


0.5072 


0.5592 


0.8041 


1.5062 


3.134 


3.0730 


1.3793 


0.5784 


0.4767 


Exact 


0.5007 


0.5072 


0.5591 


0.8038 


1.5050 


3.1314 


3.0700 


1.3800 


0.5800 


0.4702 


NGAS(2nd) 


0.5007 


0.5072 


0.5591 


0.8032 


1.5030 


3.1266 


3.0650 


1.3752 


0.5752 


0.4606 






1.5356 


1.7696 


2.7387 


5.3240 


11.193 


11.0393 


4.9984 


2.0850 


0.7713 


Caswell 






1.7695 


2.7379 


5.3216 


11.187 


11.002 


5.0900 


2.1800 


0.7703 


NGAS(2nd) 






1.7694 


2.7367 


5.3177 


11.178 


11.024 


4.9910 


2.0800 


0.7553 



Finally, we consider when the wavefunction becomes aware of the transition of the potential from a siiigle-well 
system to a double- well system in the sense that the wavefunction becomes concave from convex at y = [51]. The 
second derivative of the ground state wavefunction at y = is given by 



%£;(y) 



= i^\u^yl-ll 



which changes from negative to positive at v^y^ — 1. This equation combined with Eq. 
when a the second derivatives of the wavefunction at y = changes its sign: 



(61) 



46ll determines the moment 



2{e + \) 
3e-H 1 



0.8123, A^~ 0.418. 



(62) 



This value of critical coupling is somewhat larger than the corresponding value (0.302) of Ref. ^5] in which the approx- 
imate eigenstates of double- well oscillator were analyzed by variational approximation with 5- variational parameters. 



B. Weak coupling approximations 

In this subsection, we study the weak coupling (A^ <C 1) limit of our result. We need the weak coupling limit for 
two reasons. First, the approximate energy (j56p is too complex. To compare the result with others, it is convenient 
to consider various limiting cases. Secondly, to calculate higher order invariant operators, we need to understand 
the behaviors of the Hamiltonian ([50]) in various limits of our parameters. Explicitly, we need to know the orders 

of magnitudes of the coefficients ( — + — — ^ ) - and of the quadratic and quartic parts, respectively, of the 

\iV J e 2i2-'e^ 

Hamiltonian (|50p to construct higher order correction terms. 

Before explicitly presenting the limiting case, let us first consider the Gaussian approximation limit (yo — 0) of 

our result. The ground state wavefunction, in this case, is a Gaussian. For this wavefunction, we have the root mean 

square value of the position operator, i = = ^ ^ {VSxq), for the ground state (first excited state). Therefore, 

we have Te — 1, eo.B — Ij and Re — 3 for even parity case, and Tp = 3, eo.o = l/VS; and Rq = 5 for odd parity 
case. The value of V, is determined by the Gap equation: 

(flr f - fir - XrR = 0. 



This gap equation is the same as that of the Gaussian approximation for the state \0e)- The gap equation for the 
state |0o) is different from that of the Gaussian approximation since Rq Re- This difference makes the energy of 
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the 1** excited state different from that of the Gaussian approximation. The energies for both states are 

\h 



Eo,E/o — 



fl^E/oTE/O 



1 - 



803 



R 



E/O 



The exphcit value of Eq q is slightly larger than that of the Gaussian approximation and closer to the exact value 
of the anharmonic oscillator. This result coincides with the numerical energy eigenvalues of Ref. [3] for quantum 
anharmonic oscillators, in which a New Gaussian Approximation Scheme (NGAS) was used. 

We now consider the small coupling limit of our result for positive frequency squared (A^ ^ 1, > 0) case. In 
this case, the solution of the gap equation l[57|) becomes 



1 + 



2nlT 2 \2nlT 



XrR 



oixl 



In addition, the variational equation (j60p becomes 
d 

dyo 



XrR 



XrR 

2niT 



OiXl) 



M 

T ' 



0. 



The 0{Xr) solution to Eq. ([64]) is yo = 0. Therefore, we can series expand Eq. ([64)1 around yo = 0: 

-A 



d 
dyo 



1 



3A 9A2 



2+12 60 + ^^^0 > 



Thus, the minimizing value yc of the effective potential and become 



yc 



(3A,.)^/*, n{yc) ^1 + ^/SK + SXr + OiXl/^). 



With this the energy becomes 



2Ve 



eff 



SAt' 



3A 



0{Xl). 



From Eqs. (gZl and (gSD, we have 



e = 2yl - 2y\ 



e = l-yl + 2yt + --- . 



(63) 



(64) 



(65) 



(66) 



Therefore, we find that the correction term in is of 0{£^) ^ 0{5^) and the correction term in f2 is of 0{e) ~ 0{5). 

We next consider the A^ <C 1 limit with negative frequency squared case, < 0. In the limit j/o we have 

the exact result Q. = V2 and y^ — (\/2Ar)^^ for Ar 0. Therefore, to have series expansion around the limit, we use 
the following change of variable from to z: 



^/2Xr{l + z)' 

We have z = in the limit A^ ^ l/(\/2?/o) ^ 0. Now, the solution of the gap equation ([57|) becomes 



VL ^ no 



l-(l + ng)^^(^ + ^^ ' ^1 ■ , oo2^ 



^/2 



+ ii + nl){i + '3nl) 



XrT{l + z) 

V2 



OiXl) 



V2R 



4y4(l + z)' 



where T and R are implicitly assumed to be functions of z: 



r=l + 0(e-(-^^'-)"), ^ f-^-l) =3x/2(l + z) + ^(l + z)2 + 0(e-(v^^'-)"'; 



(67) 
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For small Ar, we essentially set T = 1. The effective energy (|59|) now becomes 



(2-»§)(l + z) (l + z)A^ 
2V2 ^ 8 



-4v^f]o(f^o + 1) + (3f^o + 14f^o - 22)(1 + z)] + 0{Xl). 



For small A^, 



2 _ V2z-^. + 0(A^). (69) 



h\uj\ dz 2 



Therefore, we have z = — — and VL = \/2 — h • • • in the small coupling limit. With this value, we have 

2v 2 2 

^ — V2 — + O(A^). The explicit dependence of e and e on Ar can be obtained from Eqs. (j48|) and ((47|) : 



£ = l-^ + 0(A?), ^=P + 0(Ay^). (70) 

We thus find that, the correction to the energy starts from the terms of O(e^) ~ 0{5^). 

With these results for the small coupling cases, we conclude that the correction to the energy starts from the order 
(5^. In addition, the corrections to the frequency and other parameters starts from the order 5. 

IV. HIGHER ORDER CONSTRUCTION OF THE CREATION AND ANNIHILATION OPERATORS 

In this section, we develop a general procedure to generate perturbative series of the invariant creation and the 
annihilation operators based on the O"* order invariant operators obtained in the previous section. First, we explicitly 
calculate the 1** order creation and annihilation operators from the 0*'' order operators. Then, we write a general 
formula which gives the n*'' order invariant creation and annihilation operators from {n — 1)*'' order ones. Some 
formulas for these calculations are given in Appendix B, C, and D. 

We first consider the structure of the Hamiltonian (l50l) which can be written in the form, 



H=^H + Eo, (71) 

where Eq is the ground state energy and the operator Ti. represents the operator part of the Hamiltonian with vanishing 
ground state expectation value, 

n = + ^T„2 + ^To^ - no; no = (0| (^T," + eh,T^ + '-^T^^ |0). (72) 

Following the principle of variational perturbation theory, we use for the parameters and yo in Eqs. (|7ip and (j72p 
their zeroth order values determined in the previous section. Note that this form ([72]) resembles Hm in Eq. (flB]) 
except for small differences in the coefficients hi. The zeroth order Hamiltonian Hm is obtained by setting hi = 1 
in 7i. In the Appendix B, we show that the coefficients hi are given, in fact, by 0{6'^) numbers by using the zeroth 
order result. The quantities 

_ ejh, 1) _ e{h, 1) 

/zi - - (73) 

are also 0{S^) numbers. In addition, it is shown that, ehi and e/i2 are of 0{S) in the Appendix B. With the 
equation ([75]) . the Hamiltonian ([7T|l can be written as the sum of Hq of and the correction terms: 

H = ^ {H„. + SAH) + Eo^ ^no, (74) 

where the perturbation term AH is 

AH=in + no- Hm)/6 = hiT^ + h2T^ ■ 
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We now show that every even power polynomials of p and ij can be written in terms of the invariant operators. 
From Eqs. ([7^ and ([TO]) we write {m + n — 1) as a function of a, d^, and H: 



a + a' 



l + \ 1 + 



2(5/i2 y V a 
2ehi 



(n + no) 



T,'-eh2{T^) 



- 2 2 



2\2 



(75) 



where ]3{t) is an 0(^°) operator, 



B(t) = 2e(H + no)-^(a-at) 



Note also that the difference between the coefficients — o? is 

Expanding Tg in series of 5, we get a well defined series expansion as a function of B, 



B 



B - 5 



-B 



0(6^); 



\B\ < 



46h2 



(76) 



(77) 



(78) 



where \B\ implies the ground state expectation value of the operator B. Since each even power polynomial of tt and y 
can be written as a sum of the multiplications of T2, T^, and Tg , they can also be written in terms of A, A\ and "H. 

We now develop the procedure to obtain the higher order invariant operators from the zeroth order ones. The 
invariant operators satisfy the LvN equation, 



^dA{p(t),x(t),t) dA ,1 , J , -, 
'^^dtt^ = + [A. H^] + S[A, AH] = 0. 

where it should be noted that Hm and AH are written as a functions of A, A' , and 7i. For example, 

2 



(79) 



AH = 



2Sh2 



^2^2 



a^ti2 



a a 

One may search for the higher order invariant operators in the form: 

A = f{a,{t),AA\-H), 



(80) 



where ais denote the coefficient functions of the polynomials of A, A^ , and 7i. 

Since A, A\ and 7i satisfy the LvN equation to 0{5'^), any function f {ai{t) ^ A, A\7i) of these operators satisfies 



^^dna.it),AA\'H) ^ ^^^/(a.(^) Ait,^) ^ 
at at 



(81) 



^d,(0 



d.f{a,{t),A,A\n) 



This equation converts the operator equation, the LvN equation, into a differential equation for the coefficients 
functions a^. 
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A. Construction of the first order invariant operators 



In this subsection, we construct the first order invariant annihilation operator from the zeroth order invariants. The 
LvN equation is solved to the first order in S and the algebra satisfied by the invariant operators are obtained. The 
1** order ground state and Hilbert space are constructed based on these first order invariant operators. 

Since the zeroth order invariant annihilation operator is A., the first order invariant annihilation operator, which 
satisfies the LvN equation ([3]), should be given by adding 0{d) correction term. Considering nontrivial accumulation 
of phase factor, we write the first order invariant annihilation operator as in the following ansatz: 



where the phase operator, ^i(i), and the 1 order correction term, Aii\, are 



at + S dt'Tiiit') 



(1)' 



hw{t) + hn{t)n, 



(82) 



(83) 



The operators in Aji) are ordered as in the normal ordering. ^'(0) is an arbitrary real constant number. A nontrivial 
operator, 7i, in the phase operator vE'i(t) is included to remove the unphysical linearly increasing term which was 
present in doing similar calculation for anharmonic oscillator in Ref. [T^ . The terms of the form A and T-IA in Ai 
are generated by the phase operator. Therefore, such terms are omitted in A(^iy 

Now, we substitute the operator Ai to the LvN equation ([3]). With the help of Eq. ([5T|) . we show that, to the first 
order in 5, the LvN equation can be solved by a simple comparison of coefScients functions with the ansatz ([5^ . To 
the first order in S, the invariant annihilation operator satisfies, 

dAi{t,m,Ht)) 



ih 

oj'l'(t) — iaOt 



dt 



.,dA 
= in—— 
dt 



hsnniiA + M(i)) + ih— 



ns 



-r!Hi(^ + M(i)) 



(84) 



where = implies "the same up to 0{S^)" . Note that the zeroth order invariant operators A and A^ in Eq. are 
given by ^ = 6oe'""*o(p(t), i(<)) and A^ = 6oe~'""*a^(p(i), x(i)) as in HH). To the zeroth order in 5, the total 
derivative dA/dt of the invariant operator vanishes. Therefore, dA/dt is an 0{S^) operator. Thus Eq. (|84p reduces 
to: 



dAn 



(1) 



dt 



By using the definition of the total derivative, the LvN equation can be written as 



' dt 



nUiA- 



1 

Sh 



ih^ + [A,H] 



(85) 



= n 



^ (2aA~[A,'H]'^ 



To the zeroth order in S, the left hand side of this equation, the total derivative 



dt 



dt 



ih 



dA dn 

by simply taking the derivatives of the coefficients, because — ~ and ~ 0. Therefore, we have 

dt dt 



dAn 



(1) 



;o,o ; 1.0 jt 
"(1) + "(i)-^ 



is computed in Appendix B: 

{^o^Tl} 



- "(1) ■ -(L)'- ■ -"{l)"-^ ""(l) 

The difference 2aA — [A, H] in the right hand side of Eq 
2aA - [A, n] 



c^'^A^ 
C(i)-^ ■ 



Sboe' 



+ Ai[ehi + eh2)T^ -TI 



IS given 



(86) 



2h2 ~2iT^ + 



On 
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where we have used the fact that both of the number ehi and e/12 are of 0{S). Using the exphcit formula for B, B^, 
and {B, a}+ in Eqs. ([76]), (|D3|l . and (|D4)) we have 



2aA-[A,'H] ih2 
~ I h[H) rr 



— + 1 1 a' + I ^ - 1 1 a- 



^t2 



(87) 



— + 1 (H + no + a)a - - 1 0^(71; + no + a) 



where 



3^2 ^ Ge^ - 

— +e{n + no) + + nof 

2 a'^ 



Therefore, with Eqs. and (|57|) . the LvN equation (ISS)) becomes 



;o,o 
^1) 



= ;^(2/iioi + 2/iii7i:i + i/i(7T;)6oe*""* 
2i L 



(88) 



(89) 



ooa a'' 



Now the LvN equation has become a simple relation between the coefficients of Eq. 
of A and HA, we get 



3{H + no + a) A - e^'^^'AHH + no + a) 

By comparing the coefficients 



6e{a + no)/i2 



hii 



6e/i2 



(90) 



After inserting Eq. ([90]) into Eq. ([M]) . we obtain A^jj by simply integrating the coefficients on the right hand side of 
the equation (|89p. Formally, we may write the formula as: 



2ai - [A, Hm] 



In summary, the creation operator and the annihilation operator to 0{5) are 

it =6ate-**^(*) =fe(aUat^))e-'*i(*), ii = 6e^*^(*)ai = 6e'*^(*)(a + a(i)), 
where the phase operator is 

6e(Sft,2 



4'i(t) = ^(O) + am 



H + no + a 



The first order correction terms in Eq. (|92p are 



"(1) — 5g~iaOt 



2a 



/i(7t;) 



2 ~t 



a^ 7Y + no + a 



2-2 



(91) 



(92) 



(93) 



(94) 



7i + no + a a H 5-a 

a"* 



2st2 



The zeroth order invariant operators satisfy a generalized deformed algebra (j26p . It is an interesting question if 
this algebraic structure is preserved to the first order. To this end, we consider the algebra satisfied by the first order 
creation and the annihilation operators, A\ and Ai. The commutator of the first order operators becomes 



452 



^2^ ^ 46^(7^ + no) _ 3e-^5h2 



1 



4(H + no)2 



a 



(95) 



The anticommutator becomes 

{Ai,A\U ^ 36^ 
462 2 



e — 



2hie5 3h2e^S 



2 7 
£5" 



in + no) + %{n + nof ^^^(n + no)^ 



a 



a" 



(96) 
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where &^ = 6q(1 + The time-dependent part of the phase factor ([93]) can be understood by calculating the 

commutator of ai and Ti., 



[ai,n] = 2a 



1 



ai, 



where the form in the square bracket of the right-hand side is remarkably the same as that appears in the phase 
factor (j93p . Thus the phase operator appears to have an imprint of the energy difference of two neighboring states in 
each parity sector. 

For these commutator and anticommutator to be those of a generalized deformed oscillator, the algebra 
must be of the following form: 

AiAl = $i(7Vi + l), ijii =$i(7Vi), 

where the structure function satisfies $i(0) = 0. Assuming the operator Ti in Eq. (|72p can be written in terms of the 
first order number operator iVi, 

niNi) = 2aNi + 5{Ei + niNi + nziVf ), 



for a = 0, 1, 2, 3 from Eqs. ((95 



we can compare the coefficients of [NiY term of AiA\ with that of A\Ai 
and ([96|) . Rather than directly computing the coefficients of (iVi)", it is convenient to compare the coefficients of the 



powers of (Ti + no) in 



with those of AiA\ by using the identity: 



H{Ni + 1) + no = 2a' + 1 + 



2 V a 



where it is assumed that 'H{Ni) is a function of Ni. The coefficients of (Ti+np)^ terms do not provide any information, 
and the coefficients of {Ti. + no)^ and (7i + no) determine ni and n2: 



ni 



12e/i2'T-o 



Uehy 



a"' 



n2 



(97) 



The zeroth order term in (7i+no) expansion is automatically satisfied with these n^'s, which is a non-trivial requirement 
to be an algebra. The energy Ei is determined by the condition = 0) = 0, which leads to a cubical equation 



262 2 



3e3(5 



e — 



4e2 2hi€S 3h2e^S 



a 



2^2 



2e' 



12/i2e^'5\ ,2 l2h2eH . 



E'o'^0, (98) 



where Eq ^ no + SEi . 

Using Eqs. ([97|) and (l98|) we get the structure function <i>i(iVi) in a series form: 

$l(iVi) ^ $o(^l) + (5(01^1 + 02^l' + 03iV3), 

7i:(7Vi)-Hno = E'o + {2a + Sni)Ni+5n2N^, 



(99) 



where cii are determined to be 



Ab''{£+~2eS), 02 



166^6 



2,2 



(100) 



Here the coefficients £+ and e are 



a 



(101) 



From $1(1) — 1, we determine the explicit value of 5^, 



46^ = 



£'+ - 215 + 4e2 ■ 



(102) 
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Notably, the structure function $i(7V'i) does not get correction term of the form since (ps = 0, even though the 
operator H gets an additional correction term of due to the 0{6) calculation. This implies that the form of the 
structure function does not change except for the 0{S) correction in coefficients. 

In summary, the creation operator, and the annihilation operator, Ai, satisfy the algebraic relation, 

ijii =$i(7Vi), iiil GS$i(iVi + l), (103) 

with the structure function (|99p. Note also that we have the desired commutation relation 

[ii,^i] = ii, [AIn,] = -A\. (104) 

From $1(0) ~ 0, the number operator annihilates the ground states: 

iVi|0£/o)i = |0£/o)i*r'(0)-0. (105) 

The n*'* excited states are generated by the formula: 

\nE/o)i = ^-=%|Ob/o)i, (106) 



where 



fc=i fe=i 



Due to Eqs. (fTU5)) and PTO|) . we have 

Ni\nE/o)i - n\nE/o)i, i{nE/o\m.E/o)i = Snm- (108) 

Eqs. (|103p . p04p . and (|105[) constitute conditions for a generalized deformed algebra which possesses a Fock space of 
eigenvectors \Oe/o)ij |l£;/o)ii ■ ■ ■ i \''T'E/o)ij • ■ ■ of the number operator iVi. 

We now explicitly calculate the P* order ground state \Oe/o)i of Eqs. (|105l) and (|106p . defined by 



^l|Oi;/o)l -0. (109) 

E/0> 



The state |0_b/o)i can be written as a hnear combination of the 0*'' order states 



\Oe/o)i = Ml ^Oe/o)o +S^^g„\nE/o)oj , 



(110) 



where Mi is the normalization constant of the state \Oe/o)i ^-nd gn are constants to be determined by the definition 
of the ground state (|109p : 

ai|0£;/o)i-0 6^gnao\nE/o)o + ai\OE/o)o (111) 

71 

= Sy^^9n+i\nE/o)o\/[n + l] +ai\OE/o)o = 0, 

n=0 

where oq is given in Eq. (|37p and the 0{S^) equation 

ao|n)o~ |n-l)ov/R, (112) 
is used. The explicit form of ai\OE/o)o is given by 



0(<5^), (113) 
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where we have used 



ieS 
2^ 



H + no 



e 



^t2 



which is derived from Eqs. ([57)l . (fTOl) . and ([75| . We have also used the fact that the zeroth order excited states are 
created by aj. The 0*'' order algebra (pS)) still holds for oq because d differs from ao only in 0((5^). From Eqs. (|llip 
and (|113p . the coefficients gi for the f** order ground state (|110p become 



51 



2^ 



(MO) 



a2 J ' 



52 



e/i2(a + np) 



53 



-2mOt 




5n>3 = 0. 



(114) 



B. General iterative formula for higher order invariant operators 

We now develop an iterative formula which gives the n*'' order creation and annihilation operators from the lower 
order operators. The procedure in obtaining a general formula is parallel to the process of finding the 1*' order 
creation and annihilation operators from the 0*^ order ones. We try to find the series solution with the following 
ansatz: 



n=0 



where A(o) = A, and the phase operator, ^'(i), and the n ' order correction term, are 

oo r,i n 

n=l •' '— " 

oo 



(115) 



(116) 



fe=0 



p,g=0 



with ^'(0) being an arbitrary real constant number, and hnk{t), 6^^'^^(t), and c^^^{t) are time dependent parameters. 

The operators are ordered as in the normal ordering. A nontrivial function of operator Ti. is included in the phase 
operator '^{t) as in the first order case. The operator 7i„ in the phase operator generates terms of the form, A and 
T-U^A. Therefore, such correction terms are omitted in ^(„). 



As we notice in Eq. ([ST|) . the total derivative 

2i dA(„) 



n dt 



2i 

n 



dt dt ih 

Y,[bl;^]{A^rn'^ + cP^;^^nnAr 
p-.q 

dA 



applies to the coefficients only: 



(117) 



To determine the parameters hnk (t) , b^,^-^ (t) , and c^^^^ (t) , we need to write 
can be accomplished by truncating the LvN equation ([3]): 



(") 



dt 



terms of A(^k) with k < n. This 



n=0 



If we truncate this equation at 0((5"), we get an iterative equation for obtaining from A(^k<7i)- 

2i d^(ri) 2i 9 ' - - 



(118) 



(119) 



where 



(n — k) ; 



(«-l) 



1 



2i dA 



(n-l) 



2(7^v4)(„_2)+A"-2) 



D 



2i dA 



(0) 



(0) 



6n dt 
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Note that the right-hand side of Eq. (|119p contains operators A(fc) with k < n. Therefore, the only remaining 
calculation to get 6(„) and C(„) is a simple comparison of the coefficients of the right-hand sides of Eqs. (|117p and 
(I119p . The comparison of the coefhcients of 7i™d determines hnk{t)- After inserting this explicit results to Eq. (|119p . 
we can integrate Eq. pi9p over time to get A^n) : 



(«) 



17 

2i 



dt' 



2(H^)(„_i)-i) 



(120) 



where it is noted that only the coefficients of the invariant operators are integrated. The calculation of A(^ij from A 
and Ti. in the previous subsection is the simplest application of this formula. 



V. SUMMARY AND DISCUSSIONS 



We have developed a new variational perturbation method which can be used for both the anharmonic and the 
double well oscillator systems simultaneously. This method is based on the observation that both the anharmonic and 
the double-well oscillator systems are parity invariant and the energy eigenstates with definite parity eigenvalues are 
approximately equidistanced for both systems. Thus the method consists of finding the creation and the annihilation 
operators which are even functions of the position and momentum operators, and raise or lower the energy eigenstates 
by two steps preserving the parity eigenvalues. To do this we expand the creation and the annihilation operators in 
an expansion parameter, 5, which is defined by the ratio of the length scales of the system and the trial wavefunctions, 
and require them to satisfy the LvN equation order by order in S. 

The zeroth order solution of the LvN equation contains two variational parameters. By minimizing the energy 
expectation value with respect to the variation of these parameters, we find the zeroth order energy eigenvalues for 
both the anharmonic and the double- well oscillators as shown in the Table 1. The errors of the numerical results 
are small enough ~ 10"^% for vast range of coupling strength, which is comparable to the 2"*^ order perturbative 
results of other Gaussian based approximations. The algebraic structure satisfied by the solutions provides a complete 
Hilbert space constructed from the variationally-determined ground state. The higher order corrections are obtained 
by applying the standard perturbation theory based on the zeroth order variational result. 

The algebra of the double-well and the anharmonic oscillator systems to the zeroth order in S is given by the 
generalized deformed oscillator with the structure function of the form, 

H ~ Eo + eiNo, 
$o(^o) - i^iNo + ^2Nl 

where A'o is the number operator to the zeroth order in 5, Eq is the ground state energy, and ei and ipi are real 
constants determined by the LvN equation. An interesting fact appears if one constructs the algebra to the first order 
in S. The Hamiltonian has the form: 

H = E'„ + e[Ni + Se'^Nl 

where Ni is the number operator to the first order in 6, and i?Q and e^ are the ground state energy and constants 
to the first order in S. The energy gets corrections proportional to the square of the number operator. On the other 
hand, the structure function to the 1** order does not get Nf correction: 

where are the constants to the 0{6). This structure function is of the same structure as that of the zeroth order 
one. This shows that the form of the structure function does not change even if we consider the first order corrections 
to the zeroth order result except for the 0{6) changes of the coefficients. It is an interesting question if this property 
continue to higher orders. 
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APPENDIX A: VARIATION OF THE EFFECTIVE POTENTIAL 



In this appendix, we give the exphcit variational calculation of Eq. ([60)1 . Varying Ve//(yo) with respect to yoi we 



get 



2 dVeffjyo) 



h\Lu\Cl dyo 



n' ( XrR 



T'- 



XrR' 



4f]3 



= 0, 



where ' denotes derivative with respect to j/o- From Eq. (|57p we have 

^, ^ sign(^^)n4M + A.4f 
- sign(.;2) + 

- - 2sign(.;^)y„^n+A.fl ^^^^ _^ ^gc'^-C) 



2Ar 



(3 + 2y2 + y2(3 + ^2)^) 



sign(tj2) 



1 



where we have used 



r 

R' 



\ / eVo ± 1 



1 ± e-v' 



(3 + 2y2 + y2(3^-2)^^_ 



1 ± e-ys 

Thus, the variational equation (|Aip becomes an equation for yc 



4yo 



±- 



1 + 



3^2 - sign(cj2) ( 1 + 2|a 



1 



(2sign(a;2)^y2 ^ ^^,^) 



2n^T ) 



-2Xr 



3^3 - sign(w2) (^1 + 17 



3^3 - sign(u;2) (i + q 
3 + 2yi+yi{3 + yi)C 

1 ± e-yo 



C I 



0. 



1 ± e-W 



(Al) 



(A2) 



(A3) 



By numerical calculation, it can be shown that this equation allows a unique positive solution for i/c for every positive 

At- . 



APPENDIX B: STRUCTURE OF THE HAMILTONIAN IN 5-EXPANSION 

In this appendix, we study the explicit behaviors of the coefficients hi in 5— expansion of the Hamiltonian (jSOp . 
which can be rewritten in the form: 

H ^ — (^2^2 + -^T^ + -Y^o) + Vixo) + Eo; (Bl) 

where Eq is the ground state energy and f2, e, e, and xq are given by the values determined by the zeroth order 
variational approximation in section III. 

The At- ^ limit corresponds to the e — > limit or e — > limit depending on the sign of uj"^ . Therefore, the analysis 
of the small limit of the variationally determined Hamiltonian, which was done in subSec. (jlll B[) . provides the 
order of magnitude of the coefficients hi with respect to e and e. Explicitly, the numbers hi and /i2 are of 0{5^). For 
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example, in the \r ^ limit of the zeroth order approximation for ground state we have, 

l--- — + 0{e^), c^2>o, A. «1, 

Xh f -l + 0(e), ^2>0, A. «1, 

Note also that the differences e{hi — 1) and e(/i2 — 1) arc of 0{S): 



e{hi — 1) 1 / Ax§ e 



-e^ + 0{e-^), Lu^>0, A^<1, 



If we define hi and /12 as 



- _ - _ S{h2 - 1) 

All - 7-3 , /l2 - TTs ' (.^-^i 



Sef^ ' 2,5 



it is clear that ehi and e/i2 are of 0((5). 



APPENDIX C: OPERATOR BASIS FOR PERTURBATIVE EXPANSION 

In this appendix, we introduce an operator basis for calculational convenience. Our basic elements are defined 
as the normalized sum of all possible terms containing m factors of tt and n factors of position operators (with the 
combination of y and (y^ — Uo))- thus a totally symmetric Hermitian operator of 0{e'^) and 0(e°), 

j=0 \ / 

For two reasons, this seems to be a natural basis with which to express operators for bi-centered systems. First, 
contains positive powers of tt and y, so it is useful for constructing generalized expansions of operators as we do in 
Eq. ([6]) in function space. Note that we can expand every nonsingular operators in terms of T," , regardless of whether 
they are symmetric or not. For example, 



where b appears as a natural expansion parameter. Secondly, T," satisfies useful set of commutation and anticommu- 
tation relations. Commutator of T," with tt and y has the effect of lowering the orders of the operator and displacing 

by yo = ^^JWY 

= -y2miT;^„i, (C2) 
[T^",^] = V2mT2"-i, [^2"+^7^] = ^/2^(2n+l)eT2" + V2^n£^2("-l), 

where the last term in the second line containing e represents the effect of bi-centered property of the system. 
Anticommutator of T," with n or y has an effect of raising the operator power and displacing by yo: 

{T^, ey}+ - V2r2"+i, {T^^+i, = ^2 {2eT^^+' + eT^") , (C3) 

{r^,7r}+ = \/2r^^p 
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When one calculates the commutator and anticommutator between higher polynomials, the following identities are 
helpful: 

[A, BC + CB] = {[A, C],B}+ + {[A, B], C}+, {A, BC + CB]+ = {{A, B}+, C}+ + [[A, C],B] (C4) 
{A,BC} = [A,B]C + B{A,C} = {A,B}C - B[A,C]. 

Similar operators as T!"j with yo for anharmonic oscillator were used in Ref. [H, and these operators have played 
a central role in studying the finite-element lattice approximation, the operator ordering, the Hahn's polynomial, and 
the analysis of exact solutions of operator differential equation. 

We list some of higher order formulae for the commutators with quadratic and quartic operators, 

[Tl:X] = 4mT^"^l\ [T^:+\T^]=A^{2n+l)eTl:_^, + A^enT'J:l-;'\ 
[T^",T2] ^ -2mzT2"+i, [T2»+\r2]^-2m»(26r2"+2^£r2"_i), 
K,Ti] 2i{n^m)eT-+2i[n/2]eT--\ 

[T^^^\n] - -4mz(2eT^n' + e«') + 2^(m)3(262T,^+32 + <5r2"_3), 
and formulae for the anticommutators with quadratic and quartic operators, 

'n-4 



1 

{T-, T^}+ = 2T:+^ - 2{mUieT-tl + sTl^^) + ^(m)4T^_4, 
where (m)„ = m(TO — 1)(to — 2) • • • (m — n + 1). All of these formulae are obtained by using \G2\ and (|C3|1 . and (|C4p 



{r;^,T2}+ = 2t;^+2 - m(m - i)er;^_ 

{T2",ri}+ = 2T,2„'VV + 2mn6T^«-i 
{r2"+i,ri}+ = AeT^+,^ + 2eT^^, + 2m{2n + l)e^T^^^,^2mn5T^J^--^,^\ 



APPENDIX D: SOME FORMULAE FOR HIGHER ORDER CALCULATION 

In this appendix, we list some formulae which are needed in evaluating higher order invariant operators. 
The operator B of Eq. ([75]) is an 0((5°) operator defined by 



^ = 2e(H + no)-^(a-at). 



With this and Eq. (jC4j) , we have 



[a,S] = 2e[a,7i] + -[d,at], 

Zl 



(Dl) 



{a,S}+ = 4e(7T; + no)a + 2e[a,7T;] + iaa^ 



-{a, at}+. 



where 



[a,7i;] = 2aa-6 



4 [a, a 



^0 



2h2 2iT^ 



4e/ii 1 _ 4«(e/ii + e/i2) 2 



II ^ z/ £H \ \(s-2e hi)T„ -e/i2(7})) 

a a 



3^2 ^ . 

+ e(7^ + no) + ^(7^ + no)2 



2 En a 







Sh2 



2(e/ii)2^ 



2(5(e2/ii 



-{'H + no,To}+ 



2q 



— T' 
2 ^0 



(D2) 
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where we have used = {Tq,TI}/2 to write in terms of A, , and H. The square of B can be written as 



= Ae^H + nof + iea{H + no, a - a+}+ + ^ ({a, a^}+ -a"- a^^) . 



(D3) 



We now Ust formulae needed to compute the commutators and anticommutators for the higher order corrections: 
[a, B] ^ Aaea - 2ia (u^s + l!!iZi±Z^^ , (D4) 



{a + a',B}+ ~ 4e (H + no + a)d + d'' {H + uq + a) + ia{a^ - d'''^), 



a To 



B' 



+ ea'^{H + no) + &e^{n + no)^ - 



aHd^ 



+2iea 



(W + no + a)a - d\n + no + a) 



dii ~ (W + 2a)a, Ha^ ■:^d\H + 2a), 



r{B,a + at}+ 



-(a^ - a^^) - [(H + no + a)a + {V. + no + a) 



{H + noX}+ 



a 
1 

1 

2^ 



2e(W + no)^ + 3e2(W + no) 1 + 



+ no + 2a)a^ + a1'^(H + no + 2a) 



+ 



4ie 



(H + no + a)^a - d\n + no + a)" 



To find the structure function to 0(5), we need the commutator and the anticommutator in terms of invariant 
operators to the same order: 



1 



[a, at] 



2Sh-2 



4e2(7i: + no) Se^ttz /. 4(7i: + no) 



1 + 



(D5) 



^(a - a^) - I (a^ + a^^) + ^ ((H + no + a)a - at(7i: + no + a)) 



-{o,oT}+ ^ +£(7{ + no) + ^(H + no)^ - ^-^(W + no) 



£2 

2 



-^(7^ + no) 1+ +^ 



4iJe2/i2 



(H + no + afd - d\H + no + af 



{n + no + 2a)a^ + at2(W + no + 2a) 



2ega 



(a -at), 



[a,wi 



2aa- 



2i/ii5 4^/^2(5 rSe^a^ 



ea^CH + no) 



+ 6e2(7i: + no)^ - ya^ +iea 3(?i: + no + a)a - at(7i: + no + a) 
The commutators between the zeroth order invariant operator and the first order correction term are 



[a, a 



hi 



4i/i2 



e(aa + 2;y^at) + 



3(1-1 + no + a)a + 2at(?i; + no + a) 



(D6) 



[o(i),at] ~ -^|£(aat + 2zy2a)- ^a2 + — [2(W + no + a)a + 3at(W + no + a)l| 
a'^ I 2 a I J J 



Summing the two, we have 

[a,at^)] + [«(!), at] 

4 

ih2 ( {c? - 2av^)e 



(D7) 



a'' 



(a - at) + ^(a^ + at^) + — + no + a)a - at(w + no + . 

2 a L 



4e^ 
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The anticommutators between the zeroth order invariant operator and the first order correction term are 
{a> «(!)}+ 



a 



h{n) + + ^ H + no + a 



a 



a 

19e2 



4e^ 



2e^ 



h e + ?io) + + 



- — a' 



h\H.) ^ 1 — \H + ?io + a) 



^at2(7T! + no + 2a) 



4i/i2 



4e^ 



2e^ 



h £ + no) + + "0 



a. 



Summing the two, we have 



{a, a|^J+ + {a(i),at}+ 



(D8) 



(D9) 



(DIO) 



ihi 
2^ 



fa -an 



4ie2/i2 



7i + no + a 



+ no + 2a) a^ + a^^{n + no + 2a) 



' a - at ( 



7i + no + a 



The commutator between the first order correction term to the annihilation operator and Ti becomes 



[a(i) , U] = a^ ( Ti + no 



Aih 



(Dll) 
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